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Interference of two independently prepared ideal Bose gases is discussed, on the basis of the idea 
of measurement-induced interference. It is known that, even if the number of atoms in each gas is 
individually fixed finite and the symmetry of the system is not broken, an interference pattern is 
observed on each single snapshot. The key role is played by the Hanbury Brown and Twiss effect, 
which leads to an oscillating pattern of the cloud of identical atoms. Then, how essential is the Bose- 
Einstein condensation to the interference? In this work, we describe two ideal Bose gases trapped 
in two separate 3D harmonic traps at a finite temperature T, using the canonical ensembles (with 
fixed numbers of atoms). We compute the full statistics of the snapshot profiles of the expanding 
and overlapping gases released from the traps. We obtain a simple formula valid for finite T, which 
shows that the average fringe spectrum (average fringe contrast) is given by the purity of each gas. 
The purity is known to be a good measure of condensation, and the formula clarifies the relevance 
of the condensation to the interference. The results for T = previously known in the literature 
can be recovered from our analysis. The fluctuation of the interference spectrum is also studied, 
and it is shown that the fluctuation is vanishingly small only below the critical temperature T c , 
meaning that interference pattern is certainly observed on every snapshot below T c . The fact that 
the number of atoms is fixed in the canonical ensemble is crucial to this vanishing fluctuation. 

PACS numbers: 03.75.Dg, 03.75.Hh, 05.30.Jp 



I. INTRODUCTION 

One of the spectacular phenomena in the physics of ul- 
tracold atomic gases is interference. When two indepen- 
dently prepared Bose-Einstein condensates (BECs) are 
released and overlap, interference fringes are observed 
between them [![. This is by no means a trivial phe- 
nomenon. Imagine two sources independently emitting 
particles toward a screen. Accumulation of the particles 
on the screen does not normally result in an interference 
pattern, since the relative phase between the two wave 
functions originating from the two independent sources is 
not well defined in general, which is crucial in Young-type 
interference experiments. 

The simplest description of the interference of inde- 
pendent BECs is based on the spontaneous symmetry 
breaking 0]: the U(l) symmetry of the system is spon- 
taneously broken upon condensation and the two gases 
individually acquire definite phases. As a result, the 
relative phase between the gases becomes well defined, 
which enables them to exhibit interference. The sym- 
metry breaking, however, would be valid only approxi- 
mately, since the actual gases in typical interference ex- 
periments consist of finite numbers of atoms. In partic- 
ular, if the number N of atoms in each gas is precisely 
fixed and its phase is completely uncertain, interference 
would not be expected between two such gases in the way 
to understand the Young-type interference. 
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Javanainen and Yoo, however, showed in their semi- 
nal paper Q that, even with such gases with fixed num- 
bers of atoms, an interference pattern can be observed 
on each snapshot photo of the overlapping gases. Notice 
that many atoms in the cloud are recorded on a photo 
at once. The indistinguishability of the identical bosonic 
atoms induces quantum correlations among them, which 
result in a nonuniform distribution of the atoms, in par- 
ticular, a sinusoidal interference pattern. The density 
profiles differ from snapshot to snapshot and the appear- 
ance of an interference pattern is not definitely certain. 
According to the numerical simulation by Javanainen and 
Yoo in [3j, however, sinusoidal patterns are very typical 
among all possible snapshot profiles and interference is al- 
most certainly observed on every snapshot (see also Q). 
The interference patterns shift randomly from snapshot 
to snapshot and the superposition of many snapshots re- 
sults in an image with no interference. This is due to the 
independence of the two gases with no phase correlation. 
One of the interference patterns and one definite relative 
phase are selected by taking a photo, i.e., by measure- 
ment, and such interference revealed on a snapshot is 
called measurement-induced interference 043 ■ 

Interference of BECs has been attracting renewed at- 
tentions these years, and a variety of interference experi- 
ments have been reported (§-Hoj. In particular, the fluc- 
tuations of the interference patterns are shown to provide 
fruitful information to probe complex many-body states 
of ultracold trapped atoms [3, 191413} , and the statistics of 
interference patterns has become an interesting subject 
to study. 

The main purpose of the present work is to clar- 
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ify the relevance of the Bosc-Einstcin condensation to 
the interference of two independent BECs, on the ba- 
sis of the idea of measurement-induced interference. As 
shown by Javanainen and Yoo, higher-order correlations 
play a crucial role for the appearance of the interference 
fringes. Then, how about the condensation? How is 
Bose-Einstein condensation essential to the interference? 
This is the question we wish to address in the present pa- 
per. Recently, Rath and Zwcrgcr have shown by a simple 
argument that the visibility of the interference is directly 
related to the condensation fraction [l3j . In this paper, 
we provide another evidence of this relationship between 
the interference and the condensation. 

We consider the following setup. We prepare two in- 
dependent ideal gases of bosonic atoms trapped in two 
spatially separated 3D harmonic traps at a finite temper- 
ature T. Each gas contains exactly N atoms separately 
and is described by a canonical ensemble with the fixed 
number of atoms. These gases are then released from 
the traps, expand in 3D free space, and overlap. Wc 
are interested in the interference patterns appearing on 
snapshot photos of the cloud of overlapping gases. 

To carry out our analysis, we compute the character- 
istic functional of the statistics of the snapshot profiles 
of the cloud, valid for the whole range of temperature 
T (across the critical temperature T c ) for a large num- 
ber of atoms N. In particular, we find that the average 
strength of the interference spectrum (Fourier spectrum 
of the density profile) over all snapshots, which is related 
to the average visibility of the interference pattern, is 
simply given by the "purity" of each gas: the larger is 
the purity of the gas, the higher is the contrast of the 
interference, and no interference is expected above the 
critical temperature T c . The purity is known to be a 
good measure of condensation [lij], because it is large 
when only a few states are macroscopically occupied and 
approaches 1 when only one state is populated. This 
shows that Bose-Einstein condensation is relevant to the 
interference. 

Furthermore, we see that the fluctuation of the inter- 
ference spectrum is vanishingly small at any temperature 
below the critical temperature T < T c , while the fluctua- 
tion abruptly changes at the critical temperature T c and 
becomes nonvanishing above T c . The interference pat- 
tern with fringe contrast depending on the purity of the 
gases is typical among all possible profiles and is certainly 
observed on every snapshot below the critical tempera- 
ture T < T c . It is shown that the canonical ensemble, 
in which the number of atoms is fixed, is crucial to the 
vanishing fluctuation. 

This paper is organized as follows. We set up tools 
to study the statistics of the snapshot profiles in Sec. 
HTl which are shown to be essentially the same as the 
ones employed in 0, [HI, [l3j| . In Sec. IIII1 the characteris- 
tic functional characterizing the canonical ensemble of a 
single gas with a fixed number of noninteracting atoms 
in a harmonic trap is given, which is the key ingredient 
in the present analysis. From this, in Sec. IIV1 we derive 



the characteristic functional of a pair of such harmonic 
clouds and compute the full statistics of the snapshot 
profiles after the two gases are released and overlap. The 
average and the covariance of the fringe spectrum are 
then analyzed in detail in Sec. [Vj obtaining the concise 
formula for the average spectrum given by the purity, 
and the fluctuation of the fringe spectrum is investigated 
as a function of the temperature. Finally, a summary of 
the work is given in Sec. I VII and some details of the cal- 
culations, concerning the derivation of the characteristic 
functional of the snapshot profiles, the treatment of the 
canonical ensemble, and the estimation of the purity and 
the other relevant quantity, are presented in Appendices 



II. STATISTICS OF SNAPSHOT PROFILES 

First of all, we setup some mathematical tools, which 
are used in the following analysis. Suppose that there 
are a large number of identical bosonic atoms and one 
takes a photo of the cloud: the positions of the N atoms 
are recorded at once on the snapshot. The probability 
of finding the ./V atoms at positions {ri, . . . ,rjy} at an 
instant t is given by 

= j^Cri) • ■ ' ft(r N )4>(r N ) ■ ■ ■ ^(n)) t , (2.1) 

where tp(r) is the field operator of the bosonic atom, 
satisfying the canonical commutation relations 



[ip(r),ipi(r')) = 5 3 (r-r'), etc. 



(2.2) 



and (■ • ■ ) t denotes the expectation value estimated in the 
state of the cloud at time t. This probability is normal- 
ized to unity as 



r N 

/ Y[d 3 nP t (N \r 1 ,...,r N ) = l. 



(2.3) 



The probability to find M atoms among N at 
{n, .. . ,rjf} is given by 

P t {M \r u ...,r M ) 

= {N N ? I)l (^(ri) ■ ■ ■ f(r^(r tf ) • • • $(n)) t , 

(2.4) 

where the normalization is such that 
J d 3 rzP t W ( ri ,...,n,...,r N ) 



P, 



(N-l) i 



(2.5) 



Given a single configuration of the ./V atoms 
{ri, . . . , rjv}, the snapshot density profile of the cloud 
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is constructed by 



1 N 

P( r ) = J^J2a{r -n), 



(2.6) 



where g(r) is a function sharply peaked around r = 
with unit volume / d 3 r g(r) = 1, characterizing the res- 
olution of the photo, and the density profile p(r) is nor- 
malized to unity: 



J d 3 rp(r) = 1. 



(2.7) 



Notice that the positions of the N atoms, {ri, . . . , rjy}, 
differ from run to run, and the density profile p(r) 
changes from snapshot to snapshot. The average pro- 
file over all possible configurations of the N atoms (over 
all snapshots) is given by 

r N 

p(r)= / Y[d 3 nP t W (r 1 ,...,r N )p(r) 

= J d 3 r'g(r-r')P t {1 \r'), (2.8) 

which gives the single-particle probability distribution 
Pf X \r), convoluted with the resolution function g(r). 

When two independent Bosc gases are overlapping, no 
interference fringes are observed in the single-particle dis- 
tribution Pj (r), which represents the image obtained 
by accumulating and superposing many snapshots. The 
result of the "independence" in the ordinary sense is that 
normal (Young-type) interference is absent, as is clear 
from (|2.8[) . An interference pattern, however, would be 
found on each snapshot, due to higher-order correlations 

We define observables that can characterize interfer- 
ence, i.e., quantities that measure whether interference 
is observed or not starting from some initial state. If 
fringes are present on a snapshot, we expect the density 
deviation 



Sp(r) = p(r) - p{r) (2.9) 
to oscillate, giving rise to spikes in its Fourier transform 



6p(k) = J d 3 rSp(r) 



(2.10) 



A spike at kf in this quantity corresponds to a spatial 
oscillation with fringe spacing 2k /kf. 

Notice here that the phase (spatial offset) of the inter- 
ference pattern varies randomly from snapshot to snap- 
shot. This is actually unavoidable, in order to be consis- 
tent with the independence of the two gases: this random 
shift smears out the fringes in the average profile p(r), 

i.e., in the single-particle distribution P ( ^(r), and the 
"independence" is recovered. In order to discard this 
random phase, we look at the square modulus of the 



Fourier spectrum, \Sp(k)\ 2 . If sinusoidal patterns with a 
definite fringe spacing (with their random spatial offsets 
discarded) are typical among all possible snapshot pro- 
files and are found on almost all snapshots, the spikes in 
the spectrum |5p(fc)| 2 would remain even in its average 
over all possible realizations of {t*i, . . . , risr}, 



S t (k) = \Sp(k)\ 2 = \p(k)\ 2 - P(k) 



(2.11) 



The typicality is characterized by the variance, or more 
generally, by the covariance 



C t (k,k') = \Sp(k)\ 2 \Sp{k')\ 2 -\Sp(k)\ 2 -\Sp(k')\ 2 . (2-12) 

If the average spectrum St(k) exhibits a nontrivial spike 
with a vanishingly small covariance Ct(k,k'), the sinu- 
soidal pattern corresponding to the spike is expected to 
be observed on every snapshot. 



By noting p(k) = g(k) £\ =1 



N -ik- 



l /N, one realizes 



that these quantities are controlled by few-particle dis- 
tribution functions. Indeed, 



|p(fc)P = |g(fc)| 2 '— 



N N 



(2.13) 



\p(k)\ 2 \p(k'W 

= \m\ 2 \g(k')\ 2 ( ^ A 7„ 1) 4 )! ^ (4) ^ fc, ) + Q (^)) 



(2.14) 



where 



l[ 2 \k) = J d 3 r 1 d 3 r 2 P t {2) (r 1 ,r 2 )e lk < ri - r2 \ (2.15) 

r 4 

4 4) (fc,fc')= / n^^ ( Vi^2,T- 3 ,r 4 ) 

a ik-(ri — r2)+ik -(rs—ri) ^ lg) 



xe 



Namely, the average fringe contrast of the N particles 

(2) 

is essentially ruled by the two-particle distribution P t , 
while its fluctuation by P^ . We do not need to compute 
the Af-particle distribution function P^ in practice to 
discuss the average fringe spectrum and the fluctuation. 
We set g{r) = 5 3 (r) henceforth, which does not spoil the 
following arguments. 

The most important feature in the two-particle prob- 
ability distribution is the Hanbury Brown and Twiss 
(HBT) effect [l5|. Due to the bosonic nature of the 
atoms, the wave function has to be symmetric under 
the exchange of particles. For instance, when two atoms 
with opposite momenta k and — k overlap, P' 2 ' oscillates 
with a period 2ir/k. In fact, in the setup discussed by 
Javanainen and Yoo fH, the initial state of the clouds 
of bosons is formed by plane waves and P { - 2 \xi,x 2 ) = 
[1 + cos k{x\ — X2)}/2. They showed that this HBT corre- 
lation is sufficient to describe the appearance of interfer- 
ence between the clouds, even if the clouds are indepen- 
dent. Inserting this two- particle probability distribution 
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in the above equations reproduces analytically their nu- 
merical result. 

In general, the fluctuation of the snapshot profiles p(r) 
is fully characterized by the generating functional 



Z t [$\ = e ' / *V *(r)p(r) _ 

When N 3> 1 , it is cast into (Appendix [5J 



(2.17) 



(2.18) 



where : • ■ ■ : denotes normal ordering. These are our tools 
for the following argument, which are essentially the same 
as the ones employed in jHHEl. 

Let us look at a Gaussian state, characterized by a 
Gaussian characteristic functional 



(W[J,J*]) t = e 



2iVNRc{a t \J) -N(J\P t \J) 



(2.19) 



where 



W[J,J*\ = e i $ d3r Jir^Hr)^ f d3r J*( r )4>(r) (2.20) 

is the normally-ordered Weyl operator, and (a t \J) = 
J d 3 ra* t (r)J(r), (J\f t \J) = f d 3 r d 3 r' J* (r)T t (r,r') x 
J(r'). In this case, the generating functional for the den- 
sity profile, Z t [$] , is given by (Appendix |A")) 



Z t [$] = e-*f dSr ^* {r) ^{W[J,J*]) i 



j,j*=o 



Det[l - i$(r) F t ] 

Tr £ (\at)(a t \+Tt/n)t<S>(r)[T t i'S>(r)] n - 1 



(2.21) 



where we have introduced an abstract notation by 

a t (r) = (r\a t ), T t (r,r') = (r\f t \r f ), (2.22) 

and 

r\r)=r\r), (2.23) 
(r\r') = 5 3 {r-r'), J d 3 r \r) (r\ = 1. (2.24) 



III. CANONICAL ENSEMBLE OF IDEAL 
BOSONIC ATOMS IN A 3D HARMONIC TRAP 

In order to discuss the interference of two indepen- 
dent ideal Bose gases released from two separate har- 
monic traps in 3D, we need to describe the state of a 
bosonic system with a fixed number of atoms. Let us 
first consider a single gas held in a trap and see how it 
is characterized by a characteristic functional (W"[J, </*])• 
The gas consists of a fixed number, N, of bosonic atoms 
and is kept in a harmonic trap at a finite temperature T. 



We assume that the 3D harmonic trap is isotropic and 
is characterized by a trapping frequency u (generalization 
to an anisotropic potential is straightforward). Let \tp n ) 
[n = (n x , n y , n z ) with n XtVtZ = 0, 1, 2, . . .] denote the en- 
ergy eigenstates of this harmonic trap. These eigenstates 
form a complete orthonormal set of bases, 



(tPn\<Pn>) = S r 



2lVJ„)< V „|=l, (3.1) 



and the field operator ip(r) is expanded as 

^(r) =^2d n if n (r), (3.2) 

n 

with d n satisfying the canonical commutation relations 

[fin,a„»] = <W, etc. (3.3) 
The Hamiltonian of the system reads 

H = E n a^a n , e„ = ^2 hum, (3.4) 



i—x,y,z 



and the number operator 



(3.5) 



When the gas is cooled below the critical temperature 
T c , the ground state \ipo) is occupied by a macroscopic 
number of atoms. In the regime 

fiw/fc B T<Cl with (huj/k B T) 3 N finite, (3.6) 

which is relevant in the actual experiments, the conden- 
sation fraction A is well approximated by (0, [l(| and 
Appendices [B] and [Cj 



,0 (T>T C ), 
with the critical temperature given by 



(3.7) 



hw f N 



k B VC(3) 



1/3 



(3.8) 



where ks is the Boltzmann constant and £(z) the Rie- 
mann zeta function. 

Since the number N of atoms in the gas is fixed, the 
gas should be described by a canonical ensemble, which 
is characterized by the characteristic functional defined 
by 



(W[J,J*}) N = 



Tr{W[J, J*]P N e~^ 6 } 
Tr{P N e~^} : 



(3.9) 
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where Pn is the projection operator onto the iV-particlc 
sector and (3 = 1/kgT the inverse temperature. In the 
regime (|3.6|) . it is estimated to be ([13] and Appendix IB"]) 



(W[J, J*]) N 



' J (2VN\(a\J)\)e- N ^ f '^ (T <T C ) 



-N(J\F\J) 



(T>T C ), 
(3.10) 



where 



)=V% ), ^ = X>«)/(e*)<P»|, (3.11) 



(T < T c ), 



\a)(a\ + P 

J={ 1„ (3.12) 



(3.13) 



TV 

with the Bosc distribution function 

The chemical potential fi (< 0) is fixed by the condition 
£/(e„) =JV (3.14) 



and is vanishingly small for T < T c . By noting a formula 
for the Bessel function 



7 ^ [ de i 
Jo 27T 



(3.15) 



the characteristic functional (|3.10[) for T < T c is equiva- 
lently expressed as 



(W[J, J*]) N 



d6 
2^ 



e 2iVNRc(a e \J) e -N(J\P'\J) 



where 



\a e ) =e l9 \a) 



(T < T c ), 
(3.16) 



(3.17) 



represents a condensate with a definite phase 6, and the 
characteristic functional (|3.10j) for T < T c is an incoher- 
ent mixture of the condensed states with different phases 
9. 

Notice that the grand canonical ensemble yields (Ap- 
pendix [B]) 

(3.18) 

with f given in p.!2p and TV being the average number 
of atoms, i.e., a different characteristic functional from 
the one for the canonical ensemble p.lO|) below the criti- 
cal temperature, while they coincide above. It is possible 
to apply the formula (|2.21[) for both canonical (|3.16[) and 



grand canonical (|3.18j) ensembles to obtain the generat- 
ing functionals for the density profile, Z t [$]. Both yield 
the same average profile [the single-particle distribution; 
see 



p(r) = T(r,r) 



(3.19) 



over the whole temperature range, while they exhibit dif- 
ferent fluctuations below the critical temperature T <T C . 



IV. TWO INDEPENDENT IDEAL BOSE GASES 
RELEASED FROM TWO SEPARATE 
HARMONIC TRAPS 

Next, we describe the release and free expansion of the 
gases and derive the state just before the measurement. 
Let us consider two ideal Bose gases independently pre- 
pared in two spatially separated harmonic traps in 3D. 
We assume that the two harmonic traps are of the same 
shape, characterized by the same trapping frequency u>, 
but spatially shifted by vectors ±d/2. The same number 
of atoms are put in the traps, TV atoms for each, at the 
same temperature T. No particle flow is present between 
the two traps. The gases are then released by turning off 
the trapping potential and expand in free space. We are 
going to discuss the interference between the overlapping 
gases. 

The energy eigenstates of the right and left traps are 
given by shifting the eigenstates \(p n ) introduced in the 
previous section, 



Iv 



(±)\ 



D T-s:P-d/2 



rr 



(4.1) 



where p is the momentum operator of an atom, which is 
the generator of the spatial shifts. We assume that the 
two traps are well separated, compared with the sizes of 
the gases, and the overlaps between the relevant eigen- 
functions of the two traps are negligible: (^n^l^i/ ) — 0. 
Under this hypothesis, they form a complete orthonormal 
set of bases for the present system, 

(v ( M> ) ) = s S s>s nn >, ££l<^ 8) ><¥i s) l = i, (4-2) 

s— ± n 

and the field operator ip(r) is expanded as 



s — ± n 



The annihilation operators al satisfy the canonical com- 
mutation relations 



^h) ■t^ri'^ ] — dss'finn'j 6tC. 



(4.4) 



The expansion of the gases in free space after the re- 
lease from the traps is easily implemented. Since the 
gases are ideal and noninteracting, the field operator 
evolves in the Heisenberg picture as 



s=± n 



(4.5) 



6 



with tp^J ( r ) = e2m v ip"n\r): we have only to replace 



(4.6) 



in any formulas, to introduce the time development. 

The two gases in the separate traps are independent 
and the state of the couple is just a product state. Each 
gas is described by the canonical ensemble with a fixed 
number N of atoms, and the characteristic functional for 
the couple is given by the product of the two character- 
istic functionals for the individual gases, each of which 
is given by (|3.10[) . or equivalcntly (|3.16[) . but shifted in 
space. Then, the generating functional for the density 
profile of the expanding and overlapping gases in free 
space after the release from the two traps is readily avail- 
able: by applying the formula (|2.2ip to the product state 
under the time evolution and by performing the integra- 
tions over the phases of the two condensates, we get 



Det[l - i<&{r)G' t ] 



z?m = < 



X Jr 



.,(+) I 



„(-)\ 



l <S>-i(r)-ig't 



I Det[l -**(*)&] 



where 



(T < T c ), 
(T > T c ), 

(4.7) 

^(-)'^ 



and 



„(±) 



, J" f (±) ', are defined by shifting rf3~TT|) - 

()3.12[) in space and time by the unitary transformations 
flUD) and gU), i.e., 



1^ 



e -*fe*e =F **- d / a |a>, (4.9) 



J- t (±) = e-i^e^P-^Te^P-^e^^, etc. (4.10) 



The grand canonical ensembles (|3.18[) . on the other 
hand, yield 



Det[l - i<f>(r)g t ] 



(4.11) 



over the whole temperature range. It coincides with the 
one for the canonical ensembles above the critical tem- 
perature T >T C , while they are different below T < T c . 



V. INTERFERENCE AND FLUCTUATION 

We are now ready to discuss the interference between 
the two gases released from the two harmonic traps, ap- 
plying the tools introduced in SecHIlto the state obtained 
in Sec. llVl 



Wc first concentrate on the generating functional (|4.7[) 
for the canonical ensembles, which is rewritten as 



Tr £ ( 

e n = lV 



X j 



=± 

oo 



^(a t (+) |i$(r)[^$(f)]"-V 



i L (-)( 



(T < T c ), 

(T > T c ). 
(5.1) 



By noting the series expansion of the Bessel function 

J ( x ) = l-l x i + l- x ± + ... , (5.2) 
the average profile is immediately obtained as 



1 



p(r)=g t (r,r) = ^fi + \r,r)+^ '(r,r)J. (5.3) 

This is just the sum of the density profiles of the two 
clouds, and no interference is observed in this quantity. 
This is the "independence" of the two gases. However, 
interference fringes are found on each snapshot. The av- 
erage spectrum defined in (|2.11l) is given in the present 
case by 



?(-), 



S t (k) = '&{&e**' p fte-*- r } 



s=± 



a\ s) \e ik -*\a { f s) )\\ 



(5.4) 



both below and above the critical temperature. Since the 
generic formula for the covariance (|2.12j) is too compli- 
cated, we just give its expressions for two extreme cases: 
at zero temperature (pure condensates), 

C t {k,k') 

= ^Re((aW|e ife -|^- ) )(aW|e- ife -H^- ) ) 

x (a^V fc '->< +) )(a^V fc '-l4 +) }) 

(T = 0), 
(5.5) 

and above the critical temperature (in the absence of 
condensates), 

C t (k,k') 

q=±k q'—^zk' 
q' = ±k' 

+ J2 \Tr{Qte lkf g t e- Wf -}\ 2 



q'=±k' 



(T > T c ). 
(5.6) 
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FIG. 1. (Color online) The time evolution of the interference 
spectrum Xt(k) in (|5,12[) between two pure condensates at 
zero temperature T = 0. The two condensates, each contain- 
ing N = 5 x 10 6 Na atoms, are released from two harmonic 
traps of trapping frequency lo = 1.6 kHz, separated by a dis- 
tance d = 30 pm. The condensates expand, overlap, and ex- 
hibit an interference pattern on each snapshot. The critical 
temperature of the gas trapped in this harmonic potential is 
estimated by (JSTHJ to be T c = 2.0 pK. 



A. At Zero Temperature T = 

Let us look at the zero-temperature case in detail. In 
this case, the average profile (|5.3|) is given by 

W) = l(\ol + \r)\* + \atHr)\ 2 ), (5.7) 

the average spectrum (|5.4[) is reduced to 

S t (fc) = i(| Xt (fc)| 2 + | Xt (-fc)| 2 ), (5.8) 

and the covariance (|5.5j) 

C t (fc,fc') = iRe(x t *(fe)x t *(-fe)Xt(-fe')Xt(fe')), (5-9) 
both given in terms of 

Xt (k) = (a[ +) \e- ik M~ ) ), (5-10) 

which is the Fourier transform of the interference term 
<Xt*{r)a t \r) between the two condensate wave func- 
tions a\ : (r). For the present harmonic traps, they read 



a\ >{r) 



nh(l + icut) 2 



3/4 



and 



(r^d/2) 2 /2h(l+iuit) 



(5.11) 



^ ^ — h[k 2 -\-(kujt—mujd/h)' 2 }/4:m,uj 



(5.12) 



The time evolution of Xt(k) is shown in Fig. [TJ for an 
experimentally realistic set of parameters. Sharp peaks 
grow in the average spectrum S t (k) given in (|5.8|) at 



k = ±kf, kj 



md 
~ftt' 



(5.13) 



The peaks become sharper and higher as time goes on. 
The covariance Ct(k,k') in (|5.9p . on the other hand, is 
vanishingly small for any (fc, k'Y since the peaks of Xt(k) 
and Xt(—k) are well separated Q. This means that there 
is no fluctuation in |<5/5(fc)| 2 around the average St(k) in 
(|5.8[) . and each single snapshot exhibits a profile 

5p(k) ~ I Xt (fc)e-*( fe ) + l Xt (-fc) e *(- fc ), (5.14) 

with an unknown phase 5 t (k). Note that there is essen- 
tially no overlap between Xt(k) and Xt(~ k), and that 
^p(r) is a real function, namely, 5p(k) = Sp*(—k). In 
addition, the phase 8t{k) should change randomly from 
snapshot to snapshot to be consistent with 5p(k) = 
[see (|5.7|) ]. By replacing St(k) — > 5/ = St(kf) since the 
spectrum is very sharp at k = ±fc/, and by performing 
the inverse Fourier transform of (|5.14[) , we get a snapshot 
profile 



p(r) 



p(r 
1 



Sp(r) 
( r ) e %/2 



„(-) 



r)e 



-iSj/2 



(5.15) 



which exhibits an interference pattern with perfect visi- 
bility, with a spatial offset Sf. Note that the visibility is 
essentially ruled by the height of the spectrum St(±fej) 
with its maximum 1/4. 

We stress that p(r) obtained in t)5 . 15[) is a snap- 
shot profile and not a quantum-mechanical average. We 
started with a fixed number N of atoms for each gas, 
and have kept the U(l) symmetry of the system during 
the whole calculation: the symmetry breaking is not as- 
sumed. In fact, the characteristic functional (|3.10p yields 



Mr)) = 0. 



(5.16) 



Without definite relative phase between the gases, in- 
terference would not be expected in the standard way 
we understand the Young-type interference. However, a 
sinusoidal pattern with a definite fringe spacing A/ = 
2ir/kf = ht/md with a good visibility is certainly ob- 
served on every snapshot, as a result of the effects of 
the higher-order correlations. Moreover, the vanishing 
covariance allows us to reconstruct the snapshot density 
profile as (|5.15[) . These are fully consistent with a series 
of the previous studies 0-0, El El . 



B. Interference and Condensation 

For a generic finite temperature T, it is possible to 
obtain asymptotic but explicit and concise formulas for 
the average spectrum St(k) and the covariance Ct(k, k') 
for large t. Let us first observe the asymptotic behavior 
of the wave function (p^}(r). By the method of steepest 
descent (stationary-phase approximation) , we get H, E3, 

m 
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V$M = / d 3 fee-^ fe2 */ 2 ™e ife -^ d / 2 )^(fe) ~ (^) 3/2 e-(^/^) 2 / 2 ^ n g (r T £)) ■ (5.17) 
The interference terms are then estimated to be 

<«ie->s, ~ (,*§)) (,±|)) e ™ 

- f d Z k' lb* fV =F —\ Co I fe' + —\ -tk'-ht(k T md/ht)/m 

~ J d 3 k' <p* n ,(k') <p n {k') e -ifc'-»t(fcTfc/)/«» ) ( 5 . 18 ) 
which is sharply peaked at k ~ ±fc/ with fey defined in (|5.13[) . and similarly, 

= /" d 3 k' if* nl {k')^ n (ky- lk '- htk / m e^ lk - {hk ' t/2m) (5.19) 



is sharply peaked at fe ~ 0. At these peaks, 



Applying these asymptotic behaviors to (|5.1[) , we get the 
average profile 



p(0) ~ 1, 



the average spectrum 



Recall here that Tr{.F 2 } is the "purity" of each gas, 
and the average fringe spectrum S t (kf) in (|5.22[) is given 
by the purity. The purity is vanishingly small Tr{J^ 2 } ~ 
(5.20) o in the absence of condensate above the critical temer- 
ature T > T c , while it becomes Tr{J" 2 } ~ 0(1) as the 
ground state is occupied by a macroscopic number of 
atoms below the critical temperature T < T c , approach- 
ing Tr{J" 2 } = 1 for pure condensation at T = 0. The 
purity is a good measure of condensation and is adopted 
for a criterion of BEC by Penrose and Onsager [l4| . The 

(5.21) 

formula for the average fringe spectrum St(kf) in (|5.22[) 
explicitly clarifies the connection between the conden- 
sation and the interference, and the importance of the 
condensation for the interference. The purity 



S t (k f )~-Tr{£ 2 }, 



5 t (0)~-Tr{^' 2 }, 



and the covariance 



(5.22) 



'C t (kf,k f ) ~ ^ (2Tr{^ 4 } + Tr{^ 2 } 2 - 3A^ 



Ct (0,0) ~ - (3Tr{J" 4 } + 2Tr{J"' 2 } 2 

C t (fc /; 0)^^Tr{J-' 4 } 

(5.23) 

for large t, where JF defined in (|3.12[) is the single-particle 
density operator of each gas, P (= T — |a)(a|) defined 
in (|3.1ip describes the thermal excitations, and A is the 
condensation fraction. The above expressions are valid 
over the whole range of temperature T. These are the 
main results of this paper. 



Tr{J" 2 } = A 2 + Tr{J-' 2 } 



(5.24) 



is different from A 2 by Tr{J"' 2 } ~ 0(1/N) [see (JU7J) in 
Appendix[C] , and therefore, the purity is essentially given 
by the square of the condensation fraction A 2 [13]. See 
Fig. [2j where the average fringe spectrum St(kf) is plot- 
ted as a function of the temperature T . 



The fluctuation of the fringe spectrum in (|5.23[) (rela- 
tive to the average), on the other hand, is estimated to 
be 



Ct(k f ,k f ) 
SKk f ) 



1 - 



A 4 - 2Tr{J"' 4 } 



(A 2 + Tr{.F 2 }) 2 
(0(l/N) (T<T C ), 
\l + 0{l/N) (T>T C ), 



(5.25) 
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FIG. 2. (Color online) The average and the fluctuation of 
the snapshot interference spectrum, SV(fc/) and Ct(kf,kf) 
given in (|5.22[) and (|5.23|) . respectively, as functions of the 
temperature of the gases T. The parameters are the same as 
in Fig.ri] The relevant quantities A, Tr{.F' 2 }, and Tr{.F' 4 } are 
numerically evaluated without resort to the continuum limit 
(IB5|I . The approximate analytical expressions ()C5[I , (1C7[) , and 
(|C9[) perfectly reproduce these numerical results. 



by noting that 

0(1), 
0(1/N), 



Tr{J^' 4 } 



0(l/N 2 ) (T<T C ), 

0(1/N 3 ) (T > T c ), 

(5.26) 



and Tr{.F' 2 } ~ 0(1/ N) for the whole temperature range 
[see ([Cl]), dUTJ), and ([U9]) in Appendix [C]. The fluctua- 
tion is vanishingly small below the critical temperature 
T < T c 3, while it is nonvanishing above T > T c . As 
shown in Fig. [3J the fluctuation abruptly changes at the 
critical temperature T c . In particular, the interference 
spectrum does not fluctuate at any temperature below 
the critical temperature T < T c , and in this range, the 
interference pattern with fringe contrast A is certainly 
observed on every snapshot. 

If the gases are described by grand canonical ensem- 
bles, instead of the canonical ensembles, the statistics of 
the snapshot profiles are given byZ t G [$] in gTTJ, and we 
end up with different conclusion from the above. In order 
to switch to the grand canonical ensembles, we have only 
to replace A and P -> f in (|Q2"|) and (|5T25j) . While 
the average fringe spectrum St{kf) remains unchanged, 
the variance Ct (fc/ , fc/ ) exhibits different fluctuation with 
the grand canonical ensembles: 

Ct(k f ,k f ) x | 2Tr{^ 4 } _ te + 0(l/N) (T<T C ), 
S t{k f ) ~ 1 + Tr{J"2}2 + (1/N) (T > T c ). 

(5.27) 

The fringe spectrum largely fluctuates below the critical 
temperature T < T c , in contrast to the vanishing fluctu- 
ation with the canonical ensembles in (|5.25p and in Fig. 
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The main difference between the canonical and grand 
canonical ensembles is the fluctuation of the total number 
of atoms. In the case of canonical ensembles, it is given 
by St{0) in (|5.22p . which is vanishingly small compared 
to the average p(0) = 1 in (|5.2ip . over the whole tem- 
perature region [Tj|. In the case of the grand canonical 
ensembles, on the other hand, it is given by 

. 1 - 2l \\\ 2 + 0{l/N) (T<T C ), 
S t (0)~-Tr{.F 2 }^ I 2 

[0(1/N) (T>T C ), 

(5.28) 

and the total number of atoms becomes fluctuating be- 
low the critical temperature T < T c . Although usually 
the canonical and grand canonical ensembles coincide in 
the thermodynamical limit N — > oo, it is not the case 
in the presence of condensate. This difference leads to 
the difference in the fluctuation of the fringe spectrum in 
(|5.25[) and (|5.27[) . From a mathematical point of view, 
the Bcssel function Jo characteristic in the generating 
functional Z^[^] for the canonical ensemble (|4.7p leads 
to the suppression of the fluctuation in the fringe spec- 
trum below the critical temperature T < T c in (15.251) . 
The canonical ensemble, in which the total number of 
atoms is fixed, is important for the interference pattern 
to be certainly observed on every snapshot below the 
critical temperature T < T c . 



VI. SUMMARY 

We have studied the interference of two independently 
prepared ideal gases of bosonic atoms, on the basis of the 
idea of measurement-induced interference. The number 
of atoms in each gas, N, is fixed finite, and the U(l) sym- 
metry of the system is not broken. Interference fringes 
are however observed on each snapshot, as a result of 
the higher-order correlations among the identical parti- 
cles. In this paper, we are interested, in particular, in 
the relevance of the Bose-Einstein condensation to this 
phenomenon [l3[. Each gas with the definite number of 
atoms N is described by a canonical ensemble trapped in 
a 3D harmonic potential at a finite temperature T [Eq. 
(|3.10p ]. and the characteristic functional -^[^ for the 
statistics of the snapshot profiles of the cloud of the over- 
lapping gases released from two spatially separated traps 
is computed [Eq. (|4.7p ], A concise formula is then ob- 
tained which clarifies the relationship between the Bose- 
Einstein condensation and the interference: the average 
fringe spectrum St(kf) is given by the purity TrjJ 72 } of 
each gas [Eq. (|5.22|) ]. which in turn is a good measure 
The fluctuation of the fringe spec- 



(E2l 

of condensation [lj| . 
trum is also analyzed [Eq. (|5.23|) ]. and the fluctuation is 
shown to be vanishingly small below the critical temper- 
ature T < T c 0] while it is nonvanishing above T > T c 
[Eq. (|5.25p and Fig. [5] . For this vanishing fluctuation, the 
canonical ensemble (the fact that the number of atoms 
in each gas is fixed) is important. 
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In the present paper, as well as in most of the previous 
works, the measurement-induced interference has been 
studied with ideal Bosc gases. It is an important subject 
to clarify the effects of the intra-atomic interaction in the 
gases . It is important to keep in mind that the single- 
particle density operator T is well defined even for an 
interacting gas, and even in this case the purity TrjJ-" 2 } 
is a good measure for the degree of condensation of the 
gas [14| . It would be tempting to apply the same reason- 
ing as the present analysis for interacting gases and to 
see how the interaction affects the fringe contrast and its 
fluctuation. However, in the case of strongly interacting 
systems, the generating functional can be substantially 
different from the one considered here. It is also an in- 
teresting problem how to deal with the interaction during 
the expansion of the gases in a self-consistent way, at least 
at the initial stages. These issues deserve investigation. 
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Appendix A: Generating Functional for Snapshot 
Profiles 

In this appendix, we sketch the derivation of the for- 
mula (|2.18p for the generating functional for the snapshot 
profiles Zt[&\ and its application to the Gaussian state 
characterized by the Gaussian characteristic functional 
(W[J, J*]) t in flUTS) to obtain Z t [$\ in (gH]) . 

The snapshot density profile p(r) in (|2.6[) fluctuates 
from snapshot to snapshot, since the configuration of the 
atoms {r 1; . . . , r^r} changes from run to run according 
to the probability distribution P t ■ The statistics of 
the snapshot profiles is characterized by the generating 
functional Z t [$] defined in (|2.17|) . i.e., 

r N 

Ztm= Hd 3 nPl N \r 1 ,...,r N )e i S d3 ^^l 
i=i 

(Al) 

Setting g(r) = S 3 (r) without loss of the essence of the 
discussion, it is arranged in the following way to obtain 
the formula in ([2~T8]l : 



■*(r iM ) (iV»l) 



(A2) 



= e -i&*dr e itf« t +c4J) e -J*r t J i (A3) 



N ■ N 

-k E *(n 

i = l 



Ztm= f f[d 3 r e P t W (r 1 ,...,r N )e 1 

oo / ■ \ M , JV 



M=0 \ / •> i=i 

M m , M 







Ml 




JL, 




MX 






ii 


M\ 





M=0 

oo ^ / . \ M - M 



(:e*r / <j3r '*( r ')V' + ('-)'/'( I ") : \ t . 



f f[d 3 r t P t {M \r 1 ,...,r M Wr 1 ) 



/ Yl d 3 r e (V3 f (ri) • • ■ ft(r M )i>(r M ) • ■ ■ 4>{ri))t 

J i 



I 

For Gaussian States is computed as 
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where the formula is written in a matrix representation 
a\j = J d 3 ra* t (r)J(r), Jt F t J = J d 3 r d 3 r' J* (r) x 
P t {r-,r')J(r'), etc., and @(r,r') = $(r)5 3 (r - r') is 
a diagonal matrix. By expressing the Gaussian fac- 
tor e~ ji:FtJ in its (path-integral) Fourier representa- 
tion, it becomes easy to apply the functional derivative 
e ~ l Tj®Tjr ^0 it, and we proceed as 



1 



Det F t 

1 / > 752 a , f; -7 t ^ 1 7 f ,i(7-at) t S , (7-at) 

Dot 1 1 



x e -i^ t (7-at) e -s ; (7-"t) t J 



Det(l - i$T t ) 

x e " 



jt(J7 1 -i*)- I J p »jt(l_iF t *)-i 0te ic4(l-i#^- t )- 1 J 

(A4) 



By putting J, J* = 0, we get (f2?2T]) . 



Appendix B: Canonical Ensemble 

The characteristic functional for the canonical ensem- 
ble (W[J, J*])n defined in (|3.9[) is to be estimated on 
the sector with a definite number TV of atoms specified 
by the projection operator Pm- It is not easy to carry 
out such a calculation in a straightforward way, but still, 
it is possible to obtain the formula for (W[J,J*])n, as 
demonstrated in [13] for the ideal Bose gas in free space. 
In this appendix, we derive the formula (|3.10l) for the 
canonical ensemble of the ideal Bose gas trapped in a 
single harmonic potential, in the regime (|3.6[) relevant in 
the ordinary experiments. 

Observe first that, by noting that 



E^ 



(Bl) 



JY 



the average in the canonical ensemble (W[J, J*])jy in 
(|3.9|) is, in general, related to that in the grand canonical 
ensemble (W[J,J*]) G in (|3~T8]) as 



N 



E 

N 

E 



Tr{W[J,J*]P N e-M"-> 1 W} 
0liN Tr{W[J,J*]P N e-P 6 } 



,/9/iAT 



Tr{ e -/3(ff-MJV)} 
Tr{P N e-^} 



N 



Tr{e-P( H -» N '>} 



(W[J,J*]) N 

(B2) 



Note that the chemical potential fi in the grand canon- 
ical ensemble (• ■ • )g is fixed by imposing 



(TV) G =]T/( £n )=TV 



(B3) 



for a given TV, where /(e) is the Bose distribution func- 
tion defined in (|3 . 1 3[) . For the ideal Bose gas trapped 
in an isotropic 3D harmonic potential characterized by a 
trapping frequency u (see Sec. IHIjl . the condition (|B3j) 
reads vM 



I 



d 3 u 



1 



(B4) 



in the limit 

fifth) -> keeping a = (/3fiw) 3 TV finite. (B5) 
There exists a solution fi (< 0) to this equation only when 

a < a = [ d 3 u — 1 

J U X ,Uy ,U Z >0 

namely, when 

T >T C (B7) 

with T c defined in ([3-8]l (but with TV instead of TV). In the 
other temperature region T < T c , the chemical potential 
is vanishing fi = and the excess number 



C(3), (B6) 



<jq = a — o = \o (T <T C 



(B8) 



is attributed to the condensed atoms, with the conden- 
sation fraction A presented in (|3.7p . 

Now, the characteristic functional for the grand canon- 
ical ensemble is estimated to be 

(W[J,J*]) G 



= e -E„/(^)|Jr>| 2 



J n = / d 3 r V * n (r)J{r) 



cxp -ao\Ja\ 



d 3 u- 



\Ju 



u x ,u v ,u z ^>0 



exp 



d 3 u- 



gUa; +Uy +U Z — \ J 

(T < T c ) 
\J U ? \ 



U X ,Uy ,U Z >0 



, — pfJ,pUx+Uy+U X 



= W G [J,J*} 
in the limit (|B5|) . where J is scaled as 

J n = ((3hLu) 3/2 J u , ii = (3hwn. 
As for the factor (Pn)g, by noting that 



-i8N 



N 



G=n[ i +/( £ ")( i 



iS M-l 



')] 



- 1 / 

(T > T c ) 
(B9) 

(BIO) 
(Bll) 
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wc have 
(phaj)- 3 (P N ) G 

K dO 



1 



-oo 2tt 1 + i£,a 



(a -a)— e-t-"-^!"* 



2tt 



,»f(<T-<r) 



<y(tr - cr) 



/C(a), 



(t < r c ) 

(B12) 



with cr = (f3fruj) 3 N . The characteristic functional for the 
canonical ensemble in the continuum limit 



(W[J,J'])n->W*[J, J*} 
is then available by inverting the relation 



(B13) 



W G [J,J*}= dalC(a)W a [J,J*}. (B14) 
Jo 

Above the critical temperature T > T c , it is just given 



by 

W a [J,J*] = W G [J,J* 



exp 



d 6 u- 

Jby 



u x ,u„,u z >0 & ° - 1 - 



(T > T c ), 
(B15) 

where ^ is fixed by (|F34|) with cr in place of cr. Below 
the critical temperature T < T c , on the other hand, the 
relation (|B14I) is essentially the Laplace transformation, 



/>oo 

/ da' e~ a ' /,Jo Wa'+a[J, J*} 
Jo 



aQe -<ro\Jo\ exp 



d A U 



u x ,u v ,u z >0 



(T < T c ), 
(B16) 

which is inverted, by noting a formula for the Besscl func- 
tion 

ds 

(B17) 



Cb 



2m 



to yield 



W„[J,J*} = J (2Vo^\Jo\ 
x exp I 



d 3 u- 



\J U 



U X ,Uy ,U Z ~>0 



oUx+Uy+Uz _ Y 



(T < T c ). 
(B18) 



Equations (|B15|) and (|B18[) arc presented in (|3.10[) , keep- 
ing in mind the limit (|B5|) and the scaling (|B10|) . 



Appendix C: Condensation Fraction and Purity 

Let us estimate Tr{.F' 2 } and Tr{J^' 4 }, which control 
the average spectrum S t (k) and the covariance C t (k, k') 
in Sec. lVBl We start by looking at Tr{P}. For the ideal 
Bose gas trapped in a single harmonic potential, setup in 
Sec. Ell we have 



N E e /3(e„-M) _ 1 



oo 



= EE< 

n j=l 



OO 



j (l_g-j/J^)3- 



(Cl) 



where z = e^ M . By splitting the occupation of the ground 
state, 

= T^ + Y z3 (t, ^W^-lV (C2) 

1-z ^ \(1 - e-^^) 3 J 



Now, since /3hu) <C 1 in the regime (13.61) , wc approximate 
it by [Hi 



Z 1 X ^ Z J 



where 



3=1 J 

Therefore, the condensation fraction is given by 
A = 1^^1-^L = 1-^V^S, (C5) 



iVl 



Tc) C(3) 



where T c is defined in (|3.8|) . This is nothing but (|3.7|) . 
Remember the condition for the chemical potential (|3.14[) 
and = C(3). 

Quite similarly, the purity of the gas is estimated as 

1 



E ( e /9(e„-/i) _ 1)2 
oo 

^ u ; fl-e- 



1 - z 



(1 - e -^ huj ) 3 

oo 

Ec?'- 1 )* 

i=i 



(1 - e -^ hul ) 3 



1 



1 - (/?^) : 



;[S2(z)-ff3 0)], 



(C6) 
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1 ' N \Tj C(3) 



(C7) 



Note that 32(1) = C(2) = n 2 /6. Furthermore, 

1 00 

£ fegg^a - 1)4 - Eo- - - 2)0- - 3)^7^ 



3=1 



z \ 1 

+ 



1-zy 6(/3ftw) 3 \i - 2 



+ 61n(l - z) + llg 2 (z) - 6g 3 (z)} 



and hence, 



Tr{ J" 4 } ~ A 4 



T 



6C(3)AT 2 \T C 



3 r 



1 + 6^ ln(l - z) J + ^[llg 2 (z) - 6.93(2)] 



(C8) 



(C9) 
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